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ABSTRACT 

Aims. To study the molecular composition of the interstellar medium (ISM) surrounding an Active Galactic Nucleus (AGN), by 
making an inventory of molecular species and their abundances. To establish a chemical differentiation between starburst galaxies and 
AGN. 

Methods. We used the IRAM-30m telescope to observe the central 1.5-2 kpc region of NGC 1068, covering the frequencies between 
86.2 GHz and 115.6 GHz. Using Boltzmann diagrams, we calculated the column densities of the detected molecules. We used a 
chemical model to reproduce the abundances found in the AGN, to determine the origin of each detected species, and to test the 
influence of UV fields, cosmic rays, and shocks on the ISM. 

Results. We identified 24 different molecular species and isotopologues, among which HC3N, SO, N 2 H + , CH3CN, NS, 13 CN, and 
HN 13 C are detected for the first time in NGC 1068. A comparison of the abundances in the nuclear regions of NGC 1068, M 82 and 
NGC 253 allowed us to establish a chemical differentiation between starburst galaxies and AGN. H2CO and CHsCCH, two abundant 
species in starburst galaxies, are not detected in NGC 1068, probably because they are destroyed by UV fields or shocks. On the other 
hand, species such as CN, SiO, HCO + and HCN, are enhanced by cosmic ray radiation fields. We obtained the upper limits to the 
isotopic ratios 12 C/ 13 C=49, 16 0/ ls O=177 and 32 S/ 34 S=5. These ratios are much lower in this AGN than in starburst galaxies. Our 
chemical models suggest that the chemistry in the nucleus of NGC 1068 is strongly influenced by cosmic rays, although high values 
of both cosmic rays and far ultraviolet (FUV) radiation fields also explain well the observations. C-shocks can explain the abundances 
of C2H and H2CO, but do not strongly affect the abundances of the other detected species. 

Conclusions. The gas in the nucleus of NGC 1068 has a different chemical composition as compared to starburst galaxies. The distinct 
physical processes dominating galaxy nuclei (e.g. C-shocks, UV fields, X-rays, cosmic rays) leave clear imprints in the chemistry of 
the gas, which allow to characterise the nucleus activity by its molecular abundances. 

Key words. ISM: molecules - galaxies:ISM- galaxies: individual: NGC 1068- galaxies: nuclei - galaxies: active 



1. Introduction 

NGC 1068 is one of the cl osest (D=14.4Mpc, 1" = 72 pc, 
Blan d-Hawthorn et"aI1 1 19971) Seyfert 2 galaxies. It has an in - 
frared (IR) luminosity of 3 x 10 u L o dTelesco & Harpedll980t) . 
IR observations reveal warm dust with temperatures of 320 K 
surrounding a smaller hot structure dJaffe et alJl2004 assumed 
to be a central black hole of mass 1.7x1 7 M ( Schin nerer et al.l 
2000). CO observations show molecular spiral arms forming a 
starburst ring 10" from the nucleus. Inside the ring, NGC 1068 
has a circumnuclear ( CND) disk with a thickness of lOpc 
(ISchinnerer et al.ll2000l see Fig. Q). 

NGC 1068 is therefore one of the best extragalactic targets 
to study the physical and chemical properties of the interstel- 
lar medium (ISM) in the vicinity of an active galactic nucleus 
(AGN). The molecular material close to the central engines of 
AGN is pervaded by X-ray and cosmic ray radiations originating 
in the nuclear accretion disks (IMalonev et al.l ll996). Previous 
single-dish and interferometric observations of NGC 1068 in a 
number of molecules, such as CO, SiO, HOC + and CN, have 



revealed a chemistry interpreted as the result of a giant X-ray 
dominated region (XDR ) in its nucleus, where shocks could also 
be heating the gas (lUsero et al.ll2004tlGarcfa-Burillo et al.ll2010T: 
Kamenetz kv et al.l 1201 1 ). This strong X-ray radiation has also 
been claimed to be responsible for the enhanced HCN-to-HCO + 
line ra tio, which differs between AGN and starburst environ- 
ments dKohno etalJl200lt iKrips et alJ l2008). However, a com- 
prehensive study of the molecular content of the ISM close to 
AGNs has not been carried out so far. As a consequence, the in- 
fluence of X-rays and cosmic rays on the molecular gas is only 
partially understood. 

Unbiased molecular line surv eys done toward th e nuclei of 
the starburst galaxie s NGC 253 dMartin et al.l 120061) and M82 
(Alad ro et al.ll201 lbh have allowed to determine the impact of 
large-scale shocks and UV fields on the ISM, as well as to 
analyse how its composition changes as the starburst phenom- 
ena evolve with time. It is clear from these studies that a de- 
tailed molecular inventory description of a prototypical AGN is 
also needed in order to understand the properties of extragalac- 
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Fig. 1: (a) 12 CO (1 - 0) map of the NGC 1068 nucleus (taken from lSchinnerer et al.ll2000|) . The beam sizes at the lowest and highest 
frequencies of the survey are plotted with a solid and dotted circles respectively. The starburst ring and CND are also indicated, (b) 
Our molecular line survey in the 3 mm wavelengths. Temperatures are in Tmb(K). 



tic XDRs, and also to chemically differentiate between nuclear 
powering sources (starbursts versus AGN). 

In this work we present a comprehensive chemical study of 
the gas in the NGC 1068 nucleus. We characterise its physical 
and chemical properties, and highlight some important differ- 
ences with the massive star-forming regions in starburst galax- 
ies. In Sects. [2] and [3] we present details of the observations, 
data reduction and analysis. Sect. [4] shows the results obtained 
from the survey, and a discussion on the NGC 1068 chemistry. 
In Sect. |5] we present a chemical model of NGC 1068 used to 
investigate how UV radiation, cosmic rays and C-type shocks 
may affect its gas. We also discuss the origin of each detected 
species and present the comparison of the model results with the 
observational data. In Sect. [6] we compare the properties of the 
circumnuclear gas in NGC 1068 with those found in starburst 
galaxies. Finally, we present our conclusions in Sect. [7] 

2. Observations and data reduction 

The observations were carried out with the IRAM-30m tele- 
scop^ (Pico Veleta Observatory, Spain) between October 2009 
and July 2010. We observed the nucleus of NGC1068, at a nomi- 
nal position Q'2ooo = 02:42:40.9, 6 2 ooo- -00:00:46.0, between the 
frequencies 86.2 and 115.6 GHz. The assumed heliocentric sys- 
temic velocity was lllOkms -1 . The half-power beam width 
(HPBW) ranged from 29" to 21", corresponding to a spatial 
scale of 2-1.5kpc. We used the band E0 of the EMIR receiver 
and the WILMA autocorrelator. This receiver-backend configu- 
ration allowed us to cover 8 GHz simultaneously in the vertical 
and horizontal polarisations, and led to a channel width spacing 
of 7 - 9 km s _1 . The IRAM-30m data were first calibrated to the 
antenna temperature (TV ) scale using the chopper-wheel method 
dPenzias & BurruslT9 73). The observations were done wobbling 
the secondary mirror with a switching frequency of 0.5 Hz and 
a beam throw of 1 10" in azimuth. We checked the pointing ac- 
curacy every hour towards the nearby bright continuum sources 
2251 + 158 and 0113-118. The pointing corrections were always 
better than 4". The focus was also checked at the beginning of 
each run and during sunsets. 

1 IRAM is supported by INSU/CNRS (France), MPG (Germany) and 
IGN (Spain). 



The observed spectra were converted from to main beam 
temperatures (Tmb) using the relation Tmb = (F e ff/B e ff)T^, 
where F s s is the forward efficiency of the telescope, whose val- 
ues were between 0.94 and 0.95, and B e ff is the main beam ef- 
ficiency, ranging from 0.77 to 0.81. Linear baselines were sub- 
tracted in all cases. We note that some of the spectra still contain 
faint ripples that were not completely removed by our switching 
scheme. This results in some of the weak features observed (e.g. 
around 94 GHz, 111.3 GHz and 1 12.3 GHz). However, these fea- 
tures are quite below the noise of the spectra. The rms achieved 
is < 2mK (< lOmJy/beam) across the whole survey. The sur- 
vey coverage is shown in Fig. |l(b)| while more detailed spectra, 
where the detected species are labeled, are shown in Fig. [2] 

The data were also corrected to first order for beam dilu- 
tion effects as T B = [(0 2 + 2 ) / 0]] J M b, where T B is the 
source averaged brightness temperature, 6 S is the source size 
and 6b is the beam size. Based on NGC 1 068 interferometric 
observations of ' 2 CO HCN and 13 CO (e.g. lHelfer &"Blit3ll99l 
ISchinnerer et al.ll2.000h we assumed an average source size of 4" 
for all the species detected in this survey. This value is consis- 
tent with other studies focused in the CND of NGC 1068 (e.g. 
| Bavetetal.ir2009bt iKrips et al.ll2008l 120111: iKamenetzkv et all 
2011). 

Gaussian profiles were fitted to all the detected lines. The 
parameters resulting from those fits are given in Table IB. ll The 
reduction of the spectra and Gaussian profile fitting were done 
using the CLASS Q and MASSA0 software packages. In those 
cases where there were two or more molecules blended, we per- 
formed a synthetic fitting (see iMartfn et al.l T20 1 Ob lAladro et al.l 
1201 lbl for details about this method) in order to calculate the 
contribution of each line to the total profile. Each case is further 
commented in Appendix lAl 

3. Analysis of the data 

The line identification was made using the CDMS ( Muller et al.l 
I200lll2005l) . JPL (lPickettetal.ll 19981) and LOVAS (lLovasll 19921) 
catalogues, and largely follows the method explained in detail 

2 CLASS http ~: / /www. iram.fr/IRAMFR /GILD AS \ 

3 MASSA http : lldamir.iem.csic.es/mediawiki - 
1 . \2.Q/index.php/MASSAJJ ser' sJAanual 
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in I Aladro et alJ (1201 lbl) for a similar line survey done toward 
the nucleus of M 82. The image band rejection of the EMIR re- 
ceiver was better than lOdB along the band, so only 13 CO(l -0) 
and CN(1 - 0) entered from the image band, at 97.9 GHz and 
94.6 GHz respectively. The latter is blended with one unidenti- 
fied feature. No strong line is expected at ~ 94.6 GHz, and the 
CN(1 - 0) line coming from the image band cannot account for 
the total emission. On the other hand, the 13 CO (1 - 0) line com- 
ing from the image band is blended with CS (2 - 1). Since there 
are no other CS lines in our survey, and the exact rejection of the 
image band at this particular frequency is not known, we can- 
not calculate the intensity of the CS transition. Unfortunately, it 
was not observed previously, so we cannot take it from the litera- 
ture. However, other CS trans itions have been ob served towards 
the nucleus of NGC 1068 bv lBavet et all ((2009b) and we include 
them in our study. 

3.1. Contribution of the starburst ring 

As shown in Fig. |l(a)| the angular resolution of the IRAM-30 m 
telescope at 3 mm wavelengths is not high enough to separate 
the CND of NGC 1068 from the surrounding starburst ring. As 
a consequence, some of the detected molecules in our survey 
might include some emission of the starburst ring picked up the 
telescope beam. 

The J - J' = 1-0 transitions of HCO + , HNC, C lg O, 
13 CO and 12 CO are deviated from a Gaussian profile. Instead, 
they show emission in the wings, and the lines are better fit- 
ted by three Gaussian components. As an example, we show in 
Fig. IA.2I the profile of the CO(l - 0) line. The gas dynamics 
in the centre of NGC 1068, and in particular in the CND, are 
quite complex. Interferometric observations have already shown 
that the profiles of CO, HCN and HCO + lines need dual, triple 
or even quadrupole Gaussian fi ts in order to better reproduce 
their profiles (Krips et al] |201 lb . This was associated to non- 
circular motions in the central ~100p c, where part of the gas 
might be blown outward due to shocks (Krips et al. 201 1), while 
the rest of the gas is accrete d toward the AGN along the jet (e.g. 
iMiiller Sanchez et alJl2009l) . 

It is not possible, with single-dish observations, to disentan- 
gle whether the deviations from the Gaussian profiles are due 
only to the complex gas dynamics or if the starburst ring emis- 
sion is also contributing to the line shapes. Furthermore, for 
other species no cited above there are other factors which make 
even more difficult to evaluate the possible starburst ring con- 
tribution. For example, some molecules such as C2H and HCN 
have hyperfine structure that intrinsically create deviations from 
Gaussian profiles. Others, such as H 13 CO + , HC3N or SiO, are 
blended in our survey or are very weak in NGC 1068, so a mul- 
tiple Gaussian fit is not feasible. Therefore, since we cannot cor- 
rect for the starburst emission, the data we show in the follow- 
ing sections include all the radiation picked up by the telescope 
beam (i.e. CND plus starburst ring). 

3.2. LTE analysis 

Assuming local thermodynamic equilibrium (LTE) and optically 
thin emission for all the detected species, we made Boltzmann 
diagrams to estimate the rotational temperatures (T Iot ) and col- 
umn densities of the molecul es (jV mo i). A detailed explan ation 
of the method can be found in Goldsm ith & Langen (11999). The 
rotation diagrams are shown in Fig. IA.3I while the resulting pa- 
rameters are listed in Table Q] To determine the rotational tem- 



peratures from the Boltzmann plots, at least two transitions of a 
given molecule are needed. In those cases where only one tran- 
sition of a molecule lies in our 3 mm survey (CO, 13 CO, C ls O, 
CS, CN, HCO + , HNC and C 34 S), we complemented our obser- 
vations with other lines, if available, in the literature (see TableQ] 
for references). For 13 CN, H 13 CN, H 13 CO + and HN 13 C there 
are no previous data published, so we simply assumed that these 
species arise from gas with the same temperature as the main 
I2 C isotopologues. We also considered the same temperature for 
HOC + and HCO + . For the rest of the cases (C 2 H, HCO, NS, 
N 2 H + ) with neither more than one transition in our survey nor 
in the literature, we fixed a r, ot = 10 + 5 K, which is an average 
value of the rotational temperatures determined in this survey. 

As mentioned before, the rotation diagrams were done as- 
suming optically thin emission of the gas. However, some 
molecules are certainly affected by optical thickness (e.g . 
CO, HCN. lJackson et al.ll993HPapadopoulos & Seaquislll 999). 
Therefore, in these cases we are underestimating the actual col- 
umn densities. The optical depths of the most intense molecules 
in the survey are further discussed in Sect. 14.31 

4. Results 

4. 1. Detected and undetected molecules 

We detected a total of 35 lines, or groups of lines, in our sur- 
vey, which correspond to 17 different molecular species and 7 
isotopologues. Some of them were detected for the first time 
in this galaxy, namely HC 3 N, SO, N 2 H + , CH 3 CN, NS, 13 CN, 
and HN 13 C. As detailed in Sect. [2] we complemented our ob- 
servations with other detections taken from the literature (see 
Table Q] for references). From the Boltzmann diagram results, 
12 CO, 13 CO, CN and CS show two temperature components, the 
warmer component coming from the higher J - J' transitions 
compiled from the literature (marked in Fig. IA. 31 with an aster- 
isk). 

We calculated column density ratios that are claimed to 
be representative of AGNs, as opposite to starburst environ- 
ments. The NGC 1068 chemistry can be summarised as fol- 
lows: high Nhcn/Nco (~ 3 x 10~ 4 ), high N C n/N HC n (~ 3 - 
4), low N H co+/Nhcn (~ 0.4), and low N HC o+/N HO c+ (~ 48). 
T hese results are co n sistent with many previous s t udies, such 
aslUsero et all d2004l) : IPerez-Beaupuits et all d2009h : lKrips eTal] 
(20081 1201 ll) . A comparison of these and other molecular ratios 
with those typically found in starburst galaxies is presented in 
Sect.[U 

We also list in Table [2] molecules that were not detected 
in our survey, but that were successfully observed in other ex- 
tragalactic sources. We show the 3<x upper limits to their col- 
umn densities, which were derived from the rms noise level of 
the corresponding spectra. Some of the non-detected species in 
NGC 1068, such as methyl acetylene (CH 3 CCH) or formalde- 
hyde (H 2 CO), are quite abundant in starburst galaxies. These two 
molecules are presumably too faint to be detected in NGC 1068, 
as other studies, also cove ring large ranges of freq uencies, did 
not detect them either (e.g. Kamenetzk v et alj|201 ll) . We discuss 
in Sect. l5.3l possible reasons for their weakness. 

4.2. Rotational temperatures 

The rotational temperatures obtained from the Boltzmann dia- 
grams (for those species with more than one detected transition) 
range between ~3 and ~51K (see Table [1). These values are 
quite below the kinetic gas temperature (7kj n ) traced by ammonia 
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Table 1 : Column densities and rotational temperatures of the ob- Table 2: Upper limits to the column densities of some undetected 



served molecules 



species. 



Molecule 



T mt f 
K 



crrr 2 



Ref 



~^CO 

13 CO 

C 18 

C 2 H 

CN 

CH 3 OH 

HCN 

CS 

HC 3 N 

HCO + 

HNC 

HCO 

SO 

NS 

HNCO 

13 CN 

NH + 

C 34 S 

H 13 CN 

SiO 

CH 3 CN 

H 13 CO + 

HN 13 C 

HOC+ 



5.0 ±0.3 
38.4 ±6.5 
4.3 ±0.2 

24.1 ±6.3 
3.3 ±0.1 
10.0 ±5.0 

8.8 ±0.1 
51.4 ±36.0 

5.6 ±0.1 

5.0 ±0.3 

7.1 ±0.4 
33.0 ± 13.1 

7.3 ± 1.4 

4.9 ±0.1 
7.6 ±0.2 
10.0 ±5.0 

22.8 ± 18.6 
10.0 ±5.0 
29.7 ±0.1 
8.8 ±0.1 
10.0 ±5.0 
3.6 ±0.5 
5.0 ± 0.3 

3.8 ±0.2 

10.2 ±0.1 

4.9 ±0.1 
7.6 ±0.2 
4.9 ±0.1 



(3.8: 
(1.2: 
(3.2: 
(7.8: 
(1.0: 
(1.0: 
(5.0: 
(5.8: 
(2.3: 
(1.4: 
(7.6: 
(7.0: 
(6.0: 
(5.3: 
(4.9: 
(2.1: 



0.3) x 
0.3) x 
0.3) x 
2.8) x 
0.1) x 
0.5) x 
0.1) x 
4.4) x 
0.1) x 
0.2) x 
2.2) x 

4.2) x 
4.0) x 
0.2) x 
0.5) x 

1.3) x 



(3.2 ± 3.3) x 



(2.8: 
(2.8: 
(2.2: 
(1.1: 
(1.6: 
(8.7: 
(6.4: 
(3.9: 
(3.1: 
(2.1: 
(1.1: 



1.4) x 
0.1) x 
0.1) x 
0.5) x 
0.6) x 
0.5) x 
0.8) x 
0.1) x 
0.1) x 
3.2) x 
0.6) x 



TO 18- 

10 18 

10 17 

10 16 

10 17 

10 16 

10 15 

10 15 

10 15 

10 15 

10 14 

10 13 

10 14 

10 14 

10 14 

10 14 

10 14 

10 14 

10 14 

10 14 

10 14 

10 14 

10 13 

10 13 

10 13 

10 13 

10 13 

10 13 



Molecule 




Molecule 






-2 

cm 




cm -2 


H 2 CO 


< 1.1 x 10 1& 


H 2 CS 


< 2.6 x 10 14 


ocs 


< 5.4 x 10 14 


C 2 S 


< 8.9 x 10 13 


CH 2 NH 


< 5.4 x 10 14 


H0C0 + 


< 7.0 x 10 13 


CH 3 CCH 


< 4.5 x 10 14 


13 cs 


< 3.6 x 10 13 


S0 2 


< 3.6 x 10 14 


NH 2 CN 


< 2.8 x 10 13 


c-C 3 H 


< 2.9 x 10 14 







See Sect.|4]for details about the calculations. 



Table 3: Molecular line ratios in NGC 1068 and comparison with 
the isotopic ratios in other galaxies. 



c 


NGC 1068 


GC 


M82 


d 


J- J' = 1-0 








li CO/ l, CO=12.6 








H 12 CO + /H 13 CO + = 17.1 








HN 12 C/HN 13 C = 23.3 


■ 12 C/ 13 C~ 20 -25 s 


12 C/ 13 C> 138 ' 




H 12 CN/H 13 CN= 16.1 








I2 CN/ 13 CN =49.1 








C 16 0/C 18 0= 176.8" 


16 O/ 18 O = 250 c 


16 0/ 18 0> 350" 


b 


C 32 S/C 34 S = 5.2* 


32 S/ 34 S~ 22 d 


32 S/ 34 S = 20' 



a: See details in Sect. 14.31 about how this value was calculated; 
b: Using data f r om th e lit erature (see Tab le ff] for re ferences) ; c: 



c-C 3 H 2 
NH 3 



7.6 



(1.7±...)x 
(1.1 ±0.2) x 



10 14 



IWilson & Roodl d!994l) ; d: iFrerking et all £l980); e 
2010), usin g C 2 H and its two C- isotopologues; /: 



Martm et al 



2009) and Aladro 1FaTl ( l20rn3) : g: IWanniedTi98(j|) . 



Martin et al 



Molecules in boldface are new detections in NGC 1068. The values 
given in this table are already corrected by beam dilution effects, 
f For those molecules with only one rotational transition available, 
we assumed T m = 10 ± 5 K. See Sect.[2]for details. 
Last c ol umn ref ers to detec ti ons t aken from the li terature : 
Israeli J2009V iMartfn et alj (120091), c iBavet etalJ J2009bl) 

rips et all i2008l) . lUsero etaP 
were not detected in this survey , but 
were compiled from the literature: f Na kaiima et alj d201 lh and 



jPerez-BeauDuits et al] d2009h 
( 120041) . c-CUH, and NH3 



were comp: 
s lAo et al] f 



20TTI) . 



(NJU), of 80 + 20 K, for the bulk of the molecular gas dAo et all 
120111) . However, we note that T mt are always lower limits to 7\j n , 
and that the conversion from rotational to kinetic temperatures 
is not simple. Also, the data we use seem to be subthermally 
excited. 

Our results agree well with previous studies where it was 
found that low -7 CO, HCN, and HCO + l ines trace co ld gas 
(Tkin° + < 20K. lPerez-Beaupuits et al.l l2009; TfP « 30 K, Israel 



2009; T, 



CO.HCN 



50 K. ISternberg et al.lll994l) . 



4.3. Isotopic line ratios 

Isotopic ratios are commonly studied through the line ratios 
between different molecular species and their isotopic substi- 
tutions. However, molecular line ratios are often affected by 
opacity effects, selective photodissociation or/and isotopic frac- 
tionation, thus giving only lower limits to the actual values 
of isotopic ratios. Fractionation processes could affect the car- 



bon i sotopic molecul a r ratios (but not the oxy gen and sulphur 
ones, IWannierl fl980l IWoods & Willacvl l2009h . and might be 
more effective in re gions with low ki netic temperatures (such 
as cold dark clouds, Sak ai et al. I l2010l) . and may be not negli- 
gible for NGC 1068 given the low temperatures traced by some 
molecules (see Sect. 14.2b . On the other hand, in the presence 
of strong UV radiation, isotope-selective photoionization can 
also take place in the external regions of molecular clouds, 
where shelf-shielding and mutual shielding between di fferent 
carbon isotopologues can occur (Glassgo ld et al] fl985l) . This 
seems to have a negligible impact on the 12 CO/ 13 CO ratio, al- 
though it may modify the C 16 0/C 18 result (IChu & Watsonl 
1983; IWilson &Ro"odlll994 . UV fields may be strong in the 
nucleus of NGC 1068 (see our modelling results in Sect. |5]l so, 
apart from our I2 CO/ I3 CO, we cannot rule out that selective 
photoionization is affecting our lower limits to the carbon and 
oxygen isotopic ratios. Finally, opacity effects clearly influence 
our derived I2 C / I3 C ratios, as they certainly affect the main iso- 
topologues of CO, HCO + and HCN in NGC 1068. On the con- 
trary, I3 C-bearing isotopic species are expected to be optically 
thin, or only moderately optically thick. 

We used the detected transitions of CO, HCN, HNC, HCO + , 
CN, and their 13 C isotopologues to calculate a lower limit to the 
I2 C / 13 C ratio in the nucleus of NGC 1068. On the other hand, 
using the observed lines of C ie O, C ls O, C 32 S, and C 34 S, we de- 
rived limits to the 16 0/ I8 and 32 S/ 34 S isotopic ratios. Our results 
are summarised in Table|3] They were calculated by dividing the 
column densities listed in Table Q] 

Our molecular line ratio 12 CO (1 - 0)/ 13 CO(l - 0) = 12.6 is 
in accordance with previous results, where moderately optically 
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thickness of t ~ 1-2 for CO (1- 0) was found dYoung & Sanders! 
ll986HPanadonoulos et al.ll996h . The HCO + (l -0)/H 13 CO + (1 - 
0),HCN(1 -0)/H 13 CN(l -0), andHNC(l - 0)/HN 13 C(l - 0) 
ratios give similar values, of 17.1, 16.1 and 23.3, respectively. 
Among all, CN is the species giving the highest lower limit to 
the 12 C/ 13 C ratio, of 49.1. 

Using the C I6 (1 -0) and C I8 (1 -0) column densities, we 
obtained 16 0/ 18 = 50.0. This low value points to the opacity 
of bot h C 16 Q(1 - 0) and C 18 Q(1 - ) lines (having the later a 
r > 1, Papadopoulos & Seaquist 1999). However, if we calculate 



this ratio using lb O/ ls O = 



x^§, and the 12 C/ 13 C = 49.1 

obtained from CN, we obtain a higher lower limit to the ls O / ls O 
ratio of 176.8. 

Unfortunately, our CS (2 - 1) line is contaminated by the 
emission of 13 CO(l - 0) coming from the image band (see 
Sect.H Therefore, we used the CS (3-2) and C 34 S ( 3-2) detec- 
tions bv lMartfn etal] d2009t) and lBavet etal] d2009bl) to calculate 
a limit to 32 S/ 34 S, which gives a result of 5.2. 



5. Modelling of the data 

To model the molecular line emissions presented in Sect. [4] we 
used the time and depth dependent chem ical model UCL_CHEM 
dViti& Williams 1 119991: IViti et alj|2004b . The molecular emis- 
sion from the AGN of NGC 1068 arises from an ensemble of 
gas components ranging in density and temperature by several 
orders of magnitude, representing various stages of the ISM, 
from diffuse to dense gas. Hence, quantitatively modelling the 
observed line emission would require some knowledge of the 
geometry as well as the dynamics of the different gas com- 
ponents. Instead chemical and photon-dominated region (PDR) 
models can be used to qualitatively determine the origin of 
each molecule and give an indication of the relative contribution 
among the different gas components. 

The UCL_CHEM was used to first simulate the formation of 
a dark core by collapsing an atomic diffuse gas from a density 
of 100 cm 3 to a final density of 10 5 cm 4 . During the collapse 
atoms and molecules freeze onto the dust grains forming icy 
mantles. Both gas and surface chemistries are self-consistently 
computed. Once the gas is in dynamical equilibrium again, 
UCL_CHEM compiles the chemical evolution of the gas and the 
dust after a burst of star formation has occurred. So, while the 
temperature during the first phase of the modelling was kept to 
10K, during the second phase it is increased to 100-350 K and 
the icy mantles are evaporated. The chemical evolution of the gas 
is then followed for 10 7 years. In both phases of the UCL_CHEM 
model the chemical network is based on more than 2345 chem- 
ical reactions involving 205 species of which 51 are surface 
species. T he gas phase react i ons were adopted from the UM IST 
data b ase dMillar et al]ll997tlLe Teuff et al]l2000tlWoodall et all 
1200% All the molecules so far detected in galactic hot cores 
(e.g. Blake et alj[l987t iMillar & Hatcheli|[T99l lHatchell et alj 
1998; Fon tani et al.ll2007l) have been included in this chemistry 
network together with all the important and familiar molecular 
species that are known from many other studies to play a sig- 
nificant role in the chemistry of interstellar gas. The surface re- 
actions included in this model are assumed to be mainly hydro- 
genation reactions, allowing chemical saturation where this is 
possible. 

The adopted values for the free parameters of the 
UCL_CHEM model are listed in Table |4] Parameters not men- 
tioned in this table have been kept to their standard (Milky Way) 
values. We used a final density of ne = 10 5 cirT 2 and a metallic - 



ity va lue ofz = 1 .056 z H , based on previous studies ofNGC 1068 
(e.g.lZaritskv et al.lll994lTacconi et1ull994lKamenetzkv et alJ 
1201 ll) . The initial elemental abundances (C/H, O/H, N/H, S/H, 
He/H, Mg/H and Si/H) used in our study are t hose correspondin g 
to extragalactic environments (as described in B avet et al.l2 008). 
scaled to the NGC 1068 metallicity. 

To sample the likely co-existing conditions that molecules 
are experiencing in the nucleus of NGC 1068, we ran five mod- 
els. The first four (a, b, c, and d; see Table|4]for details), allowed 
us to investigate the response of the chemistry to the changes 
in FUV radiation field and cosmic ray ionization rate, (. High 
cosmic ray ionization rates can be used to qualitative ly simu- 
late XDR-like environments (see e.g. lBavet e t al. 2009a). In fact, 
only few species (OH + , OH, H2O, H30 + ) show clear differences 
in abunda nces between re gions pervaded by cosmic rays and X- 
rays (Mei ierink et alj201 lb . Since none of these molecules were 
observed in our survey, for the purpose of the modelling we as- 
sume that both enhanced f and X-rays leave similar fingerprints 
in the ISM. 

Our models c and d, which are characterized by a high £ 
of 1000 times the canonical value of the Milky Way, aim to re- 
produce chemical conditions similar to the NGC 1068 nucleus. 
On the other hand, model e aims to study the possible influence 
of C-shocks, as it was recently claimed t hat they may be af- 
fecting the gas in the CND of NGC 1068 dGarcfa-Burillo et all 
I2010L iKrips et al.l 1201 lb lKamenetzkvetalJl201ll) . To do this, 
we used the initial conditions of model c and ran a version of 
UCL_C HEM coupled with the par ametric shock model devel- 
oped by Jime nez-Serra et all d2008l) . which considers a plane- 
parallel C-Shock propagating with a velocity v (we set v = 
40kms~') through the ambient m edium. De t ails of this version 
of UCL.CHEM can be found in IViti et al.l d201ll) . The influ- 
ence of shocks in the NGC 1068 nucleus is further discussed in 
Sect. 1531 

Using the results of the models, we calculated the column 
densities of the seventeen species detected in our survey 
(carbon, oxygen and sulphur isotopologues are not included) at 
a representative time of 10 5 yr (but note that chemistry has not 
necessarily reached steady state by then), using: 



N mo \ = X moX x Ay x 1.6 x 10 21 



(1) 



where X mo i is the fractional abundance of the molecule, Ay is 
the visual extinction, and 1.6 x 10 21 cirT 2 is the column den- 
sity of hydrogen at a visual extinction of one magnitude. Note 
that the formulation above simply give an 'on the spot' approx- 
imation of the column density. We set N mo \ = 10 12 cirT 2 as the 
minimum theoretical column density to consider a species de- 
tectable. Below this value we do not take into account the results 
of the models. We included, however, some species that were 
not detected in our survey (CH3CCH, H 2 CO, S0 2 , H 2 CS), but 
whose response to UV fields, cosmic rays and shocks may help 
us understand the differences in chemistry between AGNs and 
starburst galaxies. 

5. 1 . Origin of the molecules 

Disentangling the origin of extragalactic molecular emission is 
not trivial, as often the same molecules may arise from very dif- 
ferent gas conditions (e.g. Bayet et al. 2008, 2009). For the pur- 
pose of this analysis, we considered a two-components model: 
molecular emission arising from the outer layers of molecular 
clouds, which are more exposed to external radiation fields that 
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Table 4: Main parameters of the UCL_CHEM models 



Parameter 


Model a 


Model b 


Model c 


Model d 


Model e 


Temperature (phase II) T 


200/100 K 


300/100 K 


250/300 K 


350/300 K 




Visual extinction 


2& 10 mag 


2& 10 mag 


2& 10 mag 


2& 10 mag 


* 


External UV radiation intensity 


1 Habing 


1000 Habing 


1 Habing 


1000 Habing 


1 Habing 


Cosmic ray ionization rate {() 


1.3xl0~ 17 s -1 


1.3xl0~ 17 s~ l 


1.3xl0-'V 


1.3xl0- 15 s- 1 


1.3xl0- 15 s-'s 



Models a and b, characterized by a low f , aim to explore how the chemistry changes as a function of the FUV intensity. 

Models c and d, characterized by a high ( , aim to explore how the chemistry of a XDR-like environment changes as a function of FUV 

intensity. 

* In Model e the UCL_CHEM was coupled with the parametric shock model developed by Jimenez-Se rra et al .1 d2008h . and in this case the 
temperature in phase II and the visual extinction vary as the shock passes through the gas. 
1 First value corresponds to Ay = 2 mag and second to Ay = 10 mag. 



form PDRs; and emission arising from denser gas, closer to the 
centre of the clouds. These two gas components are represented 
in our model by visual extinctions of Ay = 2 mag (PDRs) and 
10 mag (dens e gas), similarly to previou s studie s (Bave t et alJ 
2008. I2009at iMeiierink & Spaansl 120051: ISpaans & MeTi erink 
2005). We note that above a certain critical visual extinction 
(~ 2.5-3 mag, depending on the radiation field) the chemistry is 
no longer affected by UV radiation and therefore any Ay above 
that could have been chosen to represent the denser, shielded gas 
component. Hence, increasing or decreasing the Ay would sim- 
ply imply a variation in the estimated column densities by the 
same factor (see Eq.Q]). 

In Table |5j we list all molecules and we classify them as 
'tracers' of a particular gas component for each model. We found 
that C2H and CN always trace PDRs and dense gas indepen- 
dently of the strength of both UV fields and CRs (i.e. accord- 
ing to all the models). On the other hand, CO, HCN, CS, SiO, 
CH 3 CN, HNC, SO, CH3CCH, H 2 CO, and S0 2 are invariably 
found in dense gas regions. We note that many molecules can 
trace PDRs as well as the denser, shielded regions, depending on 
the strength of UV fields and CRs adopted (Table 01, highlight- 
ing once more (e.g. Bayet et al. 2009) that one can not easily 
isolate PDR-tracers as templates for modelling galaxies. These 
results generally a gree with previous observational and theo- 
retical studies (e.g. iHuggins et al.lfl984t [Ste rnberg & Dalgarnol 
ll995HBavet et al.ll2008Ll2009al:lAladro et alJl201 lal) . 

5.2. Influence of UV fields and cosmic rays on molecules 

Most molecular species are sensitive to the presence of external 
fields. We checked the variations of their column densities as a 
function of UV and CRs strengths (i.e. for the models a, b, c, 
and d). A summary of the findings is shown in Table [6] While 
the column density of HOC + is enhanced by the presence of UV 
radiation, many other species are easily dissociated by UV fields 
in the external layers of the molecular clouds (such as CO, HCN, 
or NS). On the contrary, we found that methanol (CH3OH) is 
clearly destroyed by cosmic rays, while the production of several 
other species is favoured by the presence of CR photons in the 
ISM (e.g. SiO, CN, or N 2 H + ). Of particular interest is SiO, com- 
monly considered as a shock tracer, but whose enhanced abun- 
dance in NGC 1068 can also be related to hig h rates of X-rays 
(lUsero et al.ll2004t iGarc ia-Burillo et all 1201 Oh and cosmic rays 
(this work). Finally, only C 2 H shows similar column densities 
(within one order of magnitude of variation) for all the models. 

Our results agree well with those of Ba vet et all J2008), who 
also used the UCL_CHEM model to find suitable molecular trac- 
ers of hot cores under different physical conditions. In general, 



the trends of chemica l abundances with respect to f variations 
found by Bav et et alJ (|2008) are similar to what is shown in 
our Table [6] with few exceptions, being CS t he most contrasted 
species between both works. For this species Ba vet et ail (f2008) 
predicted a decrease of its abundance with increased while 
our models predict a slight increase. This may be explained by 
the different physical parameters used in both works (other ( 
strengths, initial and final densities, visual extinction and final 
temperatures). 

On the other hand, the comparison of our models with those 
that explored the PDR chemistry in a variety of extragalactic 
environments dBavet et a l. 2009a), leads to big discrepancies re- 
garding the response of the molecular abun dances as a function 
of f and FUV radiation fields. For example, Bav et et all ((2009a) 
deduced that CN and HCO + are not affected by the variations of 
4", while our results predict enhancements of those species. Also, 
iBavet et al.1 (l2009al) found that CO, HNC, HCN and H 2 CO are 
not influenced by changes of the FUV radiation field, whereas 
our results indicate that the abundances of these molecules are 
reduced by at least one order of magnitude when the FUV field 
increases. These di sagreements are mainly due to the fact that 
IBavet et al.1 (l2009al) used the UCL_PDR model where, unlike 
the UCL_CHEM code used here, the depletion of atoms and 
molecules on to grains, and the subsequent surface chemistry, 
is not included. In other words, our present modelling assumes 
that even the gas affected by UV photons has undergone suffi- 
cient processing involving gas-grain interactions. 

5.3. Comparison with the observations 

By considering a multi-component model we succeed in quali- 
tatively reproducing the observed column densities of almost all 
the detected molecules in this survey within one order of magni- 
tude, being HOC + the only exception. We consider this a good 
match, as the column densities have relevant uncertainties due to 
the assumption of LTE conditions. In general, the observations 
are better reproduced by the models which consider a high cos- 
mic ray ionization rate, but low FUV radiation field (model c), or 
high values of both parameters (model d). Once more, this result 
supports the use of high £ rates in time- and depth-dependent 
chemical models as a good approximation of XDR-like environ- 
ments. 

However, as expected, the scenario depicted by the mod- 
els reflects clear differences in the origin and timescales of the 
species, as well as highlights some degeneracies; nevertheless 
sometimes the same degeneracies can emphasise the prevalence 
of a particular energetic process. For example CO, ubiquitously 
tracing gas at a wide range of densities, is well matched by all 
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Table 5: Origin of the molecules according to each model 



Molecule 


Model a 


Model b 


Model c 


Model d 


CO 


P& D 


D 


P& D 


D 


C 2 H 


P& D 


P& D 


P& D 


P& D 


CN 


P&D 


P& D 


P&D 


P&D 


HCN 


P& D 


D 


P& D 


D 


CS 


P& D 


D 


P& D 


D 


HCO + 


- 


- 


P&D 


D 


HNC 


P& D 


D 


P&D 


D 


HCO 


P& D 


- 


P& D 


- 


N 2 H + 


- 


- 


P& D 


D 


SiO 


D 


D 


P& D 


D 


CH 3 CN 


D 


D 


P& D 


D 


NS 






P& D 


D 


hl 2 CU 


JJ 


JJ 


no, n 




HC,N 


D 








CH 3 OH 


D 








SO 


D 


D 


D 


D 


CH3CCH 


D 


D 


D 


D 


S0 2 


D 


D 


D 


D 


H 2 CS 


D 


D 






HOC + 











P stands for PDR emission, while D stands for dense gas emission 
(evaluated at visual extinctions of 2 mag and 10 mag respectively). 
P&D indicates that the molecule may arise from both molecular 
cloud regions.'-' means that the model give a column density below 
10 12 citT 2 (theoretical detectable limit), and thus the result is not 
taken into account. 



Table 6: Molecular trends with UV fields and cosmic rays 



Easily dissociated by UV fields 



Enhanced by UV fields 

Easily dissociated by cosmic rays 

Enhanced by cosmic rays 

Insensitive 



CO, CS, HCN, CH3OH, 
HCO, SO, NS, HNC, S0 2 , 
N 2 H + , SiO, CH 3 CN, H 2 CS, 
H 2 CO, HC 3 N, CH3CCH 
HOC+, 

CH3OH, HC 3 N 

CN, HCN, HCO + , HNC, SO, 

NS, N 2 H + , SiO, S0 2 , HOC + 

C 2 H 



the models, yet models with high and low radiation field, as long 
as the cosmic ray ionization rate is high, match the observations 
better. On the other hand, methanol is much better matched by a 
model where both radiation field and cosmic ray ionization rates 
are low: we interpret this 'mismatch' in physical conditions be- 
tween CO and CH3OH as an indication that the two species, 
on average, arise from different gas; methanol clearly seem to 
arise from regions where the cosmic ray ionization rate is close 
to standard, probably not representative of the average cosmic 
ray ionization rate of the galaxy. This r ate is in fact wel l known 
to be variable within our own galaxy (Dalgarno 2006), and it 
may well vary in NGC1068 too. On the other hand, the HCO + 
observed abundance is only matched if the cosmic ray ionization 
rate is high (i.e. Models c and d), independently of the FUV field 
strength. 

Unlike for cosmic rays, there are very few species that help 
us to determine the average UV field(s) strength in NGC 1068. 
Of particular interest is HCO, which seems to get closer to the 
observed value only in models where the UV field is low, while 
there is no molecule that needs a high radiation field in order to 
match the observations. While this does not exclude a high radia- 
tion field for this galaxy, especially near the nucleus, it seems un- 



likely that a rich chemistry comprising of species such as HCO 
would survive if the average radiation field were to be much 
higher than the canonical interstellar one. 

On the other hand, C2H is better reproduced in a shocked en- 
vironment, as C-type shocks enhance its abundance up to the ob- 
servational value. On the contrary, H2CO seems to be destroyed 
by C-shocks. This, as well as strong UV fields, might be the 
reason for its non-detection in our survey. Apart from C2H and 
H2CO, the presence of shocks does not greatly affect the abun- 
dances of our observed species. In particular we note that the 
SiO column density does not change if shocks are present, sug- 
gesting that its enhanced abundance is due to X-rays/cosmic rays 
gas processing instead of shocks. Hence, although they might be 
present, our models do not require the presence of strong shocks 
in the CND of NGC 1068. It is worth underlining here, however, 
that our models only include a treatment for C-type shocks: in 
these models, once the dense gas reaches its final pre-shocked 
density (i.e. at the end of Phase I), the gas and dust temperatures 
are already high enough (250K) that many exothermic reactions 
are relatively efficient. The maximum temperature reached due 
to the shock is only few hundred kelvin higher than the initial 
temperature, and moreover it does not last for long (~ 100-200 
years) before the gas starts cooling down. Hence, chemically 
there is not much difference between this scenario and the typi- 
cal 'hot core' scenario. The introduction of J-type shocks would 
significantly alter the chemistry. 

Finally, it is interesting to note that our models show different 
recycling times among molecules. This behaviour seems to be 
related to the nature of the species rather than a dependence on 
the radiation field and cosmic rays strength. For example, the 
column densities of some species, such as CH3OH and NS, are 
better matched at early times (~ 10 3 yr), since their abundances 
quickly drop at later times by at least three orders of magnitude. 
Similarly, the column densities of some important undetected 
species (e.g. CH3CCH) drop dramatically below the detectable 
levels at later times (i.e. ~ 10 5 yr). On the other hand, CO, HCO + 
and HOC + always maintain constant values. 

6. Comparison with other galactic circumnuclear 
regions 

6.1. Isotopic ratios 

For comparison purposes, Table [3] lists the isotopic line ratios 
obtained in the Galactic centre (GC) and the inner few hundred 
parsecs of the starburst galaxy M 82. As discussed in Sect. I4.3I 
our upper limits to the 12 C/ I3 C ratio in NGC 1068 are certainly 
affected by opacity effects. Also, selective photoionization and 
chemical fractionation could be altering our carbon and oxygen 
isotopic ratios in NGC 1068. These effects might partially ex- 
plain the isotopic differences among galaxies. 

The most stringent lower limit to the 12 C / 13 C ratio obtained 
in M 82 so far (> 138, lMartfn et al.ll2010l) is more than two times 
higher than our best value in NGC 1068. However, we note that 
the M 82 ratio was calculated from C2H, while for NGC 1068 we 
are using our CN results. Both molecules can be tracing different 
gas phases in the galaxies, so it is not surprising that they show 
distinct values. If we compare t he carbon isotopic ratios of both 
galaxies using CO and HCO + dHenkel & Mauersberger 1993), 
the results are similar. 

On the other hand, the difference in carbon isotopic ratios 
between galaxy nuclei may also be related to the degree of pro- 
cessed material. 12 C and 13 C are products of primary and sec- 
ondary nuclear processing in stars respectively, and thus their ra- 
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tio depends on the nucleosynthesis of stars, on stellar evolution, 
and therefore on the star formation history of the galaxy. The 
stellar population a ge in the centre of NGC 1068 (200-300 Myr, 
Davies et aD2007h is much older than that of M 82 (~ 15 Myr, 
Konstantopoulos et al. 2009). Primary-to-secondary isotopic ra- 
tios are predicted to decrease with time, so a lower I2 C / I3 C ra- 
tio (as well as le O/ ls O) is expected in NGC 1068 than in M 82. 
On the other hand, both the Galactic centre and the NGC 1068 
nucleus have low carbon isotopic ratios as compared to M82 
(Table [3]). This implies that both nuclei have a high chemical 
processing degree. 

Our limit to the 32 S/ 34 S in NGC 1068 is a factor of four 
lower than the ratios in the GC and M 82. In the Milky Way, 
it was observed that 32 S / 34 S does not strongly vary with time, 
although increases with the galacticentric distance, following 
the 12 C/ 13 C gradient of the galactic di sk dFrerking et al.l l 1980; 
I Wilson & Rood! 1 1994 IChin et al.lfl996h . Also, it was proposed 
that this ratio depends on the stellar population, so that environ- 
ments with more massive stars and supernova explosions (which 
is the case of NGC 1068 as compared to the GC and M 82) have 
lower 32 S/ 34 S. 

6.2. Molecular abundances 

M82 and NGC 253 are the starburst galaxies with the best 
known chemical compositions so far. Like NGC 1068, both have 
been t he targets of unbiased molecular line sur veys toward their 
nuclei dMartm et al.ll2006t IXIadro et alj|201 lbl) . The starburst in 
the centre of M82 is old, with an average stellar population age 
of ~10Myr dKonstantopoulos et ai1l2009l) . This, together with 
a high supernova rate, creates strong UV fields, so that PDRs 
dominate its chemistry. On the other hand, NGC 2 53 is claimed 
to be in an earlier stage of starburst evolution ( Marti n et al.l 
120091; lAladro et alj|201 lab, with younger stellar p opulations in 
its nucleus (~ 6 Mvr. iFernandez-Ontiveros et al.l 12009b where 
the PDRs, although present, do not characterise its ISM. Instead, 
the nucleus of NGC 253 show that large-scale shocks are the 
dominant process driving the molecular abundances. 

Figure [3] shows the abundances of all species detected by 
single-dish mm/sub-mm telescopes in the three galaxies, as well 
as other species that were identified in both starburst galaxies 
but remain undetected in NGC 1068 (for them we used the upper 
limits to the column densities listed in Table|2]). This figure aims 
to highlight the chemical differences between AGN-dominated 
and starburst-dominated nuclei. The fractional abundances were 
calculated with respect to C 34 S, which is significantly less af- 
fected by opacity effects than CS, and whose abundance does 
not se em to be signific antly affected by the type of nuclear activ- 
ity dMartm et al.ll2009b . 

We note that while M 82 and NGC 253 are located at about 
the same distance (3.6 and 3.9 Mpc respectively, Freedman et al.l 
Il994t iKarachentsev et al1l2003h . NGC 1068 lies four times fur- 
ther. However, while the comparison of the column densities in 
Fig. [3] would be biased by the different distances to the galax- 
ies, this should not affect the comparison of the fractional abun- 
dances. 

Among all the species compared, only the ratios 
H 13 CN/C 34 S and H 13 CO + /C 34 S are higher in NGC 1068 
than in the starburst galaxies, being between two and five 
times more abundant in the AGN. Conversely, the fractional 
abundances of C 2 H, CH 3 OH, HC 3 N, HCO and SO show similar 
values in both types of nuclei, with differences within a factor 
of two among the galaxies. This suggests that these species 
are not especially affected by the different nuclear activities 



in starburst ga l axies and AGN, as previously suggested by 
Nakaiima et alJ (1201 ll) . A third group is formed by those species 
whose fractional abundances are clearly higher in both M 82 and 
NGC 253 than in NGC 1068. Examples are CH 3 CN/C 34 S, or 
HOC + /C 34 S, which show differences in the abundances between 
three and five times. It is remarkable that these molecules 
have very similar abundances in M82 and NGC 253. However, 
the biggest contrasts are given by those molecules that have 
been clearly detected in both starburst galaxies, but not in our 
NGC 1068 survey. Examples are CH3CCH, H 2 CO, c-C 3 H 2 , 
H 2 CS, S0 2 and NH 2 CN. We note that propyne (CH3CCH) and 
formaldehyde (H 2 CO) are abundant specie s in starburst galaxies 
dMiihle et alJ l2007t lAladro et al.l 1201 lal) . Unlike H 2 CO, we 
obtained a stringent upper limit to the CH3CCH column density 
in NGC 1068 (Table |2|, which makes it the most contrasted 
species of our survey when comparing starburst galaxies and 
AGN (Fig. [3J. We found that CH 3 CCH/C 34 S is -40 and ~5 
times more abundant in M 82 and NGC 253, respectively, than 
in NGC 1068. This strong variation in the CH3CCH fractional 
abundances is likely to be intrinsically related to the different 
nuclear activities in galaxies, rather than being due to other 
factors, such as different opacities (due to its K-ladder structure 
the transitions of this molecule are expected to be optically 
thin), or kinetic temperature of the gas (as we do not find a 
strong ( anti-) correlation of its abundances with the 7\i n of the 
three galaxies). 

7. Conclusions 

We used the IRAM-30m telescope to carry out a spectral line 
survey (from 86.2 GHz to 1 15.6 GHz) towards the nucleus of the 
AGN galaxy NGC 1068. We detected a total of 35 spectral lines 
corresponding to 17 molecular species. Seven molecules, or their 
13 C isotopologues, are detected for the first time in this galaxy, 
namely HC 3 N, SO, N 2 H\ CH3CN, NS, 13 CN, and HN 13 C. We 
note that some molecules that are abundant in other galaxies, 
such as H 2 CO or CH3CCH, are not detected in NGC 1068. We 
calculated the column densities of the identified species assum- 
ing local thermodynamic equilibrium and optically thin emis- 
sion. However, some intense lines such as J - J' = 1-0 CO, 
HCO + , HCN or HNC, are affected by optical thickness (and 
might be also altered by chemical fractionation and isotope- 
selective photoionization effects), so their column densities are 
lower limits to the actual values. 

Using several lines from our survey we obtained the isotopic 
abundance ratios 12 C/ 13 C= 49, 16 0/ 18 0> 177 and 32 S/ 34 S = 
5 for NGC 1068. These isotopic ratios significantly differ from 
those of the Galactic centre and starburst galaxies. Apart from 
the fact that the opacities of the transitions change from source 
to source, this reflects the different star formation histories of the 
galaxies. 

We ran chemical models to determine the origin of the 
species detected in our survey, as well as to test the influence 
of UV fields, cosmic rays, and shocks on their abundances. We 
used a high rate of cosmic ray ionization rate to qualitatively 
simulate a XDR-like environment. We found that the assumption 
of some molecules being PDR tracers or dense gas tracers should 
be done carefully, and only if the physical conditions of the gas 
are known, since under certain circumstances a PDR tracer may 
be also arising from dense gas regions and viceversa. Our mod- 
els indicate that a high cosmic ray ionization rate is present in 
the nucleus of NGC 1068, albeit a combination of a high cosmic 
ray ionization rate and UV fiels may also explain the observed 
abundances in the galaxy. On the contrary, shocks do not seem 
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to be necessary to explain most of the molecular abundances. 
The chemical models applied to NGC 1068 predict that while 
UV radiation easily dissociate a large number of molecules (e.g. 
HCN, CH 3 OH, HNC, or CH3CCH), cosmic rays enhance the 
abundances of a good number of others, such as CN, SiO, N2H + , 
NS, or HCN. 

We aimed to find the molecules that are more sensitive to 
the different power sources of galaxies nuclei (starbursts and ac- 
tive galactic nuclei). We carried out a detailed comparison of 
the chemical composition of the ISM in the NGC 1068 nucleus 
and those of the two starbursts galaxies M82 and NGC 253. 
Examples of differing ratios are H 13 CN/C 34 S, which is ~3-5 
times more abundant in NGC 1068 than in starburst galaxies, and 
CH3CCH/C 34 S, which is, at least, ~5-40 times more abundant in 
starburst galaxies than in NGC 1068. 
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Fig. 2: NGC 1068 frequency survey coverage between 86.2 and 1 15.6 GHz. The spectra have been smoothed to a velocity resolution 
of~44kms~ 1 (corresponding to ~ 16MHz).The assumed systemic velocity is lllOkms -1 . 
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Fig. 3: Comparison of the fractional abundances among the two starbust galaxies M82 (blue) and NGC 253 (red), and the AGN 
NGC 1068 (green). Y-axis is logio(N mo i/N C 34 S ). Molecules marked with an * were not detected in our NGC 1068 survey; arrows 
indicate their upper limits to the fractional abundances. 
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Appendix A: Comments on individual molecules 

We have detected seventeen molecular species, and seven car- 
bon, sulphur, and oxygen isotopic substitutions in the 3 mm at- 
mospheric window, between 86.2 GHz and 115.6 GHz. Here, 
we detail each observed transitions, the Gaussian profile fitting, 
cases of blended lines and other special cases. In order to identify 
the molecules, we started from the premise that only the transi- 
tions having E\ ovi < 100 cm -1 and log( J Tmb dv) > -6 nm 2 MHz 
contribute to the observed lines. After checking this initial con- 
dition, we ensured that the spectral features that not fulfil these 
conditions can be considered negligible. 

NGC 1068 has a heliocentric systemic velocity of v sys = 
1137kms _1 (value taken from NASA/IPAC Extragalactic 
Database - NED). Based on the values obtained from the 
Gaussian fits, we have chosen an average value of 1 100 km s -1 
for fixing this parameter when needed. We found two different 
line widths for the detected molecules in NGC 1068. In gen- 
eral, strong lines have values in the range FWHM~ [230 - 
260] km s _1 , while faint lines have line widths of about ~ [120 - 
190]kms _I . The cases where it was necessary to fix this pa- 
rameter are commented below, since the value depends on each 
species. 

- Hydrogen cyanide - HCN and H I3 CN 

We detected the H 13 CN(1 - 0) transition at 86.3 GHz. It 
is formed by three hyperfine spectral features, although the 
structure of the line is not resolved due to the broad line 
widths of NGC 1068. To take into account the contribu- 
tion of each hyperfine feature we have performed a syn- 
thetic Gaussian fit (details about this m ethod can be seen in 
iMartfn et al.ll2010t I Aladro et al.ll201 lbl) . The resulting val- 
ues, shown in Table IB.fl correspond to the main component 
of the group. This HCN carbon isotopologue is detected for 
the first time in this galaxy. 

HCN(1 - 0) is identified at 88.6 GHz. In a similar way as 
for H I3 CN, we fitted a synthetic Gaussian profile leaving the 
linewidth and line position free. 

- Hydrogen isocyanide - HNC and HN 13 C 

One strong line of hydrogen isocyanide was detected at 
90.7 GHz. The telescope beam possibly picked up some 
emission coming from the starburst ring that surrounds the 
CND of NGC 1068. We complement ed this transition with 
the H NC (4 - 3) line observed by |Perez- Beaupuit s~et al.l 
(200§. 

The HN 13 C(1 -0) transition at 87.1 GHz has been tentatively 
detected for the first time in this galaxy. Because it is quite 
faint, the integrated intensity resulting from the Gaussian fit 
is highly affected by the baseline. This is reflected in the col- 
umn density derived from the rotational diagram, which has 
a high associated error. 

- Ethynyl - C 2 H 

C 2 H (1 - 0), formed by six hyperfine features, is detected at 
87.3 GHz. We performed a synthetic Gaussian fit to calculate 
the contribution of each feature. Because we only have one 
detected transition of this species, we calculated its column 
density fixing the rotational temperature to 10 ± 5 K. This is 
the fourth most abundant molecule after CO and its carbon 
and oxygen isotopologues. 

- Isocyanic acid - HNCO 

Four groups of transitions were detected at 87.9 GHz, 
109.5 GHz, 109.9 GHz and 110.3 GHz. They are composed 
by 5 features each. The line at 109.9 GHz is seen as a bump 
at the higher frequencies of C ls O, although both can be 



easily separated. The transition at 110.3 GHz is very close 
to 13 CO and blended with CH 3 CN(6 4>0 - 5 4j0 ). We per- 
formed a synthetic Gaussian fit to all these lines, includ- 
ing other two HNCO transitions that are too faint to be de- 
tected, but that fulfil the conditions mentioned at the begin- 
ning of this appendix (£(„„ <100cm _I and log( J Tmb dv) > 
-6 nm 2 MHz). In the Gaussian fitting we fixed the linewidth 
to 230 km s _1 and the line position to 1 100km s . The inte- 
grated area of our HNC O (5q,s - 4p 4) line i s consistent with 
the values obtained by Marti n~et~al] (|2P09). However, from 
the Boltzmann diagram we obtain a rotational temperature 
of r rot ~ 3 K, which is almost 5 times higher than the one 
obtained bv lMartfn et al.l d2009l) . This difference is due to the 
transitions with higher energy levels we included. 

- Oxomethylium - HCO + and H 13 CO + 

The strong line of HCO + (1 - 0) at 89.2 GHz is the only de- 
tected transition in our line survey, and shows some emis- 
sion coming from the starburst ring. We als o used the 
HCO+(3 - 2) line observed at 267.6 GHz bv iKrips et all 
(2008). The H 13 CO + (l - 0) line is detected at 86.8 GHz, al- 
though it is blended with HCO(li, - U) ) and SiO(2 - 1). 
To disentangle the contribution of each molecule, we per- 
formed synthetic Gaussian fits to the three lines. We first fit- 
ted H 13 CO + (l - 0) fixing the position to 1 100 km s -1 and 
the linewidth to 240kms~ I . Then, we fitted the other two 
molecules (see Fig. lA.U . 

- Oxomethyl - HCO 

Only one molecular transition of this species has been de- 
tected along the survey. It is blended with H 13 CO + (1-0) and 
SiO (2-1) (see Figs. |3]and |ATj- HCO (1 1, - I>0 ) is located 
at 86.7 GHz, and contains four features. After subtracting the 
oxomethylium Gaussian fit, we over-fitted a Gaussian profile 
to the HCO features with a fixed position of 1 100 km s _I and 
a linewidth of 240 km s . 

- Silicon monoxide - SiO 

As commented before, SiO (2 - 1) at 86.8 GHz is blended 
with H 13 CO + (l - 0) and HCO(li, - 0i, ). After subtract- 
ing the Gaussian fits of the other two molecules, we fitted 
another synthetic Gaussian profile to the residuals, fixing 
only the position to 1100 km s~' (see Fig. IA. lb . This is the 
strongest molecule of the three blended lines. 

- Methyl cyanide - CH 3 CN 

This symmetric top molecule is detected for the first time in 
NGC 1068. The two CH 3 CN (5* - 4 k ) and CH 3 CN (6 k - 5 k ) 
transitions at 92.0 GHz and 110.4 GHz lie inside our fre- 
quency coverage. Each of them are formed by 25 spec- 
tral features with the quantum number k going from to 
/ - 1 . Although faint, both lines are strong enough to be de- 
tected with a S/N > 3. The transition at 1 10.4 GHz is blended 
with HNCO(5i j4 - 4i 3 ). We fitted a synthetic Gaussian 
profile to both lines, fixing the line position of CH 3 CN to 
1 100km s" 1 and the linewidth to 150km s _] (based on the 
non-blended line). Its rotational temperature is similar in 
NGC 1068 (T ml ~ 29 K) and M 82 (T mt ~ 33 K). 

- Sulphur monoxide - SO 

Three transitions of sulphur monoxide were detected for 
the first time in NGC 1068 at 99.3 GHz, 100.0 GHz, and 
109.252 GHz. The first one was not blended and its Gaussian 
parameters have been used to fit the other two. The line 
is narrow, 136kms~'. The other two lines are blended 
with HC 3 N(11 - 10) and HC 3 N(12 - 11) repectively. The 
contribution of SO was firstly subtracted and then we fit- 
ted cyanoacetylene to the residuals. The SO(5 4 - 4 4 ) at 
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Fig. A.l: Gaussian fits to the blended lines of H 13 CO + (1-0), 
HCO (1 1 o-Oi,o) and SiO (2-1). We first fitted a Gaussian profile 
to the h' 3 CO + line, and then removed its emission. After that, 
we fitted the hyperfine structure of HCO, and also removed it. 
Finally, we fitted a SiO Gaussian profile to the residuals. See 
more details in the text. 




Fig. A. 2: Line profile of the CO (1 -0) line. Y-xis is in T m \, scale. 
The contribution of the starburst ring emission was fitted by two 
Gaussians at the wings of the central Gaussian, which contains 
the CND emission. 



100.0 GHz does not fit well in the rotational diagram, so the 
values we show in Table Q~|do not take into account this tran- 
sition. 

Hydroxymethylidynium - HOC + 

We detected the HOC + (l - 0) transition at 89.5 GHz. Its 
peak intensity temperature is slightly above the noise, so 
we consider it as a tentative detection. This line shows 
a quite narrow linewidth, of about 128kms _1 . Fixing 
the rotational temperature to 10 + 5K, we obtained the 
lowest column density of the molecular survey, with 
AW =(l.l±0.6)xl0 13 cnr 2 . 
Cyanoacetylene - HC3N 

This molecule is detected for the first time in NGC 1068. 
We identified HC 3 N(10 - 9) at 91.0GHz, HC 3 N(11 - 10) 
at 100.1 GHz, and HC 3 N(12 - 11) at 109.2GHz. The tran- 
sition at 9 1.0 GHz was not blended, so we used its Gaussian 
parameters to fix the other two, which were blended with SO. 
For those lines, we firstly subtracted the sulphur monoxide 
features (as explained before), and then fitted two synthetic 
Gaussian profiles to the residuals. The HC3N rotational tem- 
perature in NGC 1068 (7.3 + 1.4K) is significantly lower than 
the ones found fo r M 82 and NGC 25 3 (-25 K and ~[12- 
73] K respectively. lAladro et aD201 lal) . 



Diazenylium - N2H + 

This is the first time that this species is detected in 
NGC 1068. Only the N 2 H + (1 - 0) transition, at 93.2 GHz, 
lies inside our survey. It is formed by three features, so we 
performed a synthetic Gaussian fit in order to calculate the 
contribution of each one. We fixed the rotational temperature 
to 10 + 5 K to calculate its column density. 
Carbon monosulfide - CS and C 34 S 

We detected the CS (2 - 1) line, which is blended with 
13 CO(l - 0) coming from the image band. Since these two 
are the only transitions of these molecules detected in our 
survey, it is not possible to dissentangle the contribution of 
each one. Thus, no Gaussian fit was done. Nevertheless, we 
us ed the CS (3 -2), CS (5 - 4), and CS (7 - 6) lines observed 
by iBavet et al.l d2009bl) to do the rotation diagram. Two gas 
components with different temperatures and column densi- 
ties are seen. 

C 34 S(2 - 1) was detected at 96.4 GHz. This transition 
was complemented with the C 34 S (3 - 2) line observed by 
iMartm et alj ((2009). The rotational temperature is one of the 
lowest obtained in the survey, with only ~ 4 K. 
Carbon monoxide - CO, C 18 0, and 13 CO 
CO ( 1 - 0), at 1 1 5 .3 GHz, was blended with NS (3 1 - 2 , ), but 
since NS was comparatively very faint, we assumed its con- 
tribution negligible. As the emission of CO(l - 0) is quite 
spread in the nucleus of NGC 1068, it clearly shows a devi- 
ation from a Gaussian profile, which points to the complex 
gas dynamics in that region, as well as a possible contribu- 
tion of the starburst ring emission. We fitted a three Gaussian 
profile, which is shown in Fig. IA.2l 

We a lso used the C O (2-1), CO (3 - 2), and CO (4 - 3) lines 
from llsrael d2009l) . Like CS, CO shows two components in 
the rotation diagram. Despite CO (1 - 0) is the strongest fea- 
ture in the survey, it is not seen in the image band. This is 
due to the variable rejection along the 3 mm band. 
13 CO(l - 0) was detected at 110.2GHz. Similarly to 
12 CO(l - 0), this transition showed a deviation from a 
Gaussian profile. It was compleme nted with the 13 CO (2 - 1) 
and 13 CO (3 - 2) lines observed bv llsraell d2009l) . Two differ- 
ent components are also seen in the Boltzmann diagram. As 
commented before, 13 CO (1-0) does appear in the image 
band, blended with CS (2 - 1 ). 

C ls O (1 - 0) appears at 109.8 GHz. It was not blended, but it 
showed a bump at the higher frequencies due to the proxim- 
ity of a HNCO transition. To avoid the isocyanic acid con- 
tamination, we fixed the linewidth to 230 km s 1 and then fit- 
ted three Gaussians, which better account for th e line profile. 
We als o used the C ls O (2-1) line observed by Marti n et all 
(2009) for complementing our observations. The rotational 
temperature is the lowest one among all the molecules de- 
tected (T mt = 3.3 + 0.1 K). 
Cyanogen - CN 

Two separated lines of cyanogen were identified at 
1 13.2GHz (CN (l ,i -0 ,i)) and 1 13.5 GHz (CN (1 ,2 -0 ,i))- 
The first one is formed by four different features, while the 
second one has five. We fitted a synthetic Gaussian profile 
to both lines. We included ot her CN transitions observed by 
iPerez-Bea upuits et al. (2009). We also detected two multi- 
feature lines of 13 CN at 108.4 GHz and 108.6 GHz. 
Methanol - CH 3 OH 

Only the CH 3 OH (2 k - l k ) line is seen at 96.7 GHz, which is 
composed by four features with the quantum number k vary- 
ing between — 1,0, 1. However, we calculated the integrated 
intensities of other transitions that lie inside our frequency 
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range, which fulfil the conditions imposed to the energy and 
logarithm of the integrated intensities, as explained at the be- 
ginning of this appendix. Taking into account all these fea- 
tures (thirteen in total), we performed a synthetic Gaussian 
fit. 

- Nitrogen monosulfide - NS 

We detected the NS (3i - 2i) transition in our survey. It is 
splitted in two separated multi-lines with six features each. 
The first one is blended with CO (1 -0), while the second one 
is located at higher frequencies, so that can be easily sepa- 
rated from CO. We firstly subtracted the carbon monoxide 
contribution, and then fitted a synthetic Gaussian profile to 
all the NS features. We assumed r rot = 10 + 5 K to calculate 
the column density. 
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Fig. A. 3: Boltzmann diagrams of the molecules observed in the survey for which we have detected more than one transition, or for 
which we have taken other transitions from the literature (marked with *). The resulting rotational temperatures are indicated. 
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Appendix B: Gaussian fits parameter results. 



Line 


Frequency 


/ 


Position 


Width 


Tpeak 


Comments 




MHz 


Kkms" 1 


kms~' 


km s -1 


mK 




H 13 CN(1 -0) 


86340.2 


0.4 ±... 


1066.3 + 0.0 


166.4 ±0.0 


2.4 


hf, s 


HCO(li, -0i, ) 


86670.8 


0.2 ±... 


1100.0 ±0.0 


240.0 + 0.0 


0.6 


b, s 


H 13 CO + (l -0) 


86754.3 


0.5 ± ... 


1100.0 ±0.0 


240.0 + 0.0 


1.9 


b, s 


SiO(2- 1) 


86847.0 


0.6 ±... 


1100.0 ±0.0 


259.1 ±0.0 


2.2 


b, s 


HN 13 C(1 - 0) 


87090.8 


0.2 ±0.3 


1075.7 + 124.4 


151.8 + 235.9 


1.1 




C 2 H(1 -0) 


87316.9 


3.1 + ... 


1098.4 ± 0.0 


240.0 + 0.0 


12.0 


hf, s 


HNCO(4 04 - 3 03 ) 


87925.2 


0.1 + ... 


1100.0 ±0.0 


230.0 + 0.0 


0.4 


m, s 


HCN(1 -0) 


88631.8 


11.1 + ... 


1120.1 ±0.0 


238.5 + 0.0 


43.0 


hf, s 


HCO + (l -0) 


89188.6 


13.8 + 0.3 


1122.6 ± 1.2 


230.9 ± 1.5 


55.9 




HOC + (l - 0) 


89487.4 


0.2 ±0.1 


1173.3 + 17.5 


128.2 ±51.6 


1.8 




HNC(1 - 0) 


90663.6 


7.5 ± 0.6 


1141.3 ±0.3 


235.9 + 3.2 


30.8 




HC 3 N(10-9) 


90979.0 


1.1 ±0.1 


1142.6 ± 10.0 


263.7 ± 24.9 


3.8 




CH 3 CN(5 k -4 k ) 


91987.0 


0.1 ±... 


1100.0 ±0.0 


150.0 ±0.0 


0.6 


m, s 


N 2 H + (1 -0) 


93173.7 


1.0 ±... 


1164.1 ±0.0 


233.5 + 0.0 


4.0 


m, s 


C 34 S (2 - 1) 


96412.9 


0.7 ±0.1 


1091.2 + 15.7 


249.0 ± 29.5 


2.7 




CH 3 OH(2 t - \ k ) 


96741.4 


0.9 + .. . 


1117.5 ±0.0 


247.2 + 0.0 


3.3 


m, s 


CS(2- 1) 


97981.0 










b 


SO(2! -3 2 ) 


99299.9 


0.4 ±0.1 


1111.1 ±7.5 


136.0 ± 17.0 


3.1 




SO(5 4 -4 4 ) 


100029.6 


0.2 ±0.1 


1111.1 ±0.0 


136.0 ± 0.0 


1.6 


b 


HC 3 N(11 - 10) 


100076.4 


0.5 + ... 


1142.6 ±0.0 


263.0 + 0.0 


2.0 


b, s 


13 CN(l li0 -0 M ) 


108426.9 


0.1 ±... 


1104. ±0.0 


257.2 ± 0.0 


0.2 


m, s 


13 CN(l 2 ,i -0i' ) 


108651.3 


0.2 + ... 


1104. ±0.0 


257.2 ± 0.0 


0.6 


m, s 


HC 3 N(12- 11) 


109173.6 


0.7 ± ... 


1142.6 ±0.0 


263.0 ± 0.0 


2.6 


b, s 


SO(3 2 -2i) 


109252.2 


0.2 ±0.1 


1111.1 ±0.0 


136.0 + 0.0 


1.3 


b 


HNCO(5!, 5 -4 M ) 


109496.0 


0.04 ± ... 


1100.0 ±0.0 


230.0 + 0.0 


0.2 


m, s 


C 18 0(1 - 0) 


109782.2 


3.6 ±0.2 


1142.2 ±6.1 


230.0 + 0.0 


13.0 




HNCO(5 ,5-4o, 4 ) 


109905.8 


0.3 + .. . 


1100.0 ±0.0 


230.0 + 0.0 


1.3 


m, s 


13 CO(l -0) 


110201.4 


12.8 ±0.7 


1105.2 ±3.6 


252.1 ±4.9 


45.8 




HNCO(5i4-4 13 ) 


110298.1 


0.2 + ... 


1100.0 ±0.0 


230.0 ± 0.0 


0.7 


m, b, s 


CH 3 CN(6,-5,) 


110383.5 


0.1 + ... 


1100.0 ±0.0 


150.0 ±0.0 


0.8 


m, b, s 


CN(l ,2-Oo,i) 


113191.3 


4.5 + ... 


1152.8 ±0.0 


247.4 + 0.0 


17.3 


s 


CN(l ,2-0 ,i) 


113491.0 


12.1 + ... 


1152.8 ±0.0 


247.4 ± 0.0 


46.0 


s 


NS(3!-20 


115153.9 


0.4 + ... 


1100.0 ±0.0 


140.0 ± 0.0 


2.8 


m, b, s 


12 CO(l -0) 


115271.2 


172.2 ±4.5 


1168.3 ±0.3 


251.5 ±0.4 


618.0 


b 


NS (3 1 -2 1 ) 


115556.2 


0.4 + ... 


1100.0 ±0.0 


140.0 + 0.0 


2.8 


m, s 



Table B.l: Parameters fixed during the Gaussian fit have zero errors associated. Synthetic fits do not give errors to the integrated 
area. Remarks: (b) blended line; (m) multi-transition line; (hf) hyperfine transition; (s) synthetic Gaussian fit using MASSA. For 
those transitions showing several components, the parameters refer to the main component of the group. 



2 



